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SUMMARY 


This  paper  can  be  divided  into  essentially  two  different  parts.  The 
first  part  contains  a  qualitative  discussion  of  the  one -electron  approxima¬ 
tion  and  a  rough  review  of  the  B ohm -Pines,  van  Hove -Hugenholtz  and 
Luttinger  work.  These  works  can  be  considered  as  representative  of  that 
part  of  current  many -body  literature  which  aims  directly  at  the  behavior 
of  quasi -particles  in  normal  many-fermion  systems. 

The  second  part  of  this  paper  tries  to  elaborate  a  specific  method 
to  handle  quasi -particles.  The  method  is  formally  similar  to  the  usual 
Hartree-Fock  scheme  and  it  was  chosen  out  of  many  possible  approaches 
partly  because  the  author  has  been  actively  working  with  it  and  partly 
because  it  is  fairly  simple  and  still  illustrates  the  ideas  and  scope  of  present 
day  many -body  theories. 

Some  of  the  material  given  here  has  not  been  published  earlier  to 
the  knowledge  of  the  author.  That  goes  for  the  iteration  scheme  developed 
to  approximate  the  self-energy  operator  M  ,  and  also  for  the  evaluation  and 
discussion  of  the  properties  of  the  first  non-trivial  iteration  step  M^^^  , 
for  some  different  physical  systems,  particularly  for  the  hydrogen-like 
levels  of  alkali  atoms.  It  is  further  shown  that  the  Bohm -Pines  treatment 
of  single -particle  properties  corresponds  to  an  "adiabatic'^  approximation 
of  M^^^  .  This  approximation  of  M^^^  is  however  able  to  deal  naturally 
also  with  the  periodic  lattice  case. 
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INTRODUCTION 

The  name  one -electron  approximation  has  almost  become  synonymous 
with  the  Hartree-Fock  method  (HF).  Here  we  will,  however,  use  that  name 
for  any  scheme  where  a  one --particle  equation  generates  orbitals  and  orbital 
energies  which  are  then  used  to  describe  the  properties  of  the  system. 

When  the  first  one-electron  approximation  was  presented  by  Hartree^^ 

in  1928,  it  was  received  with  considerable  reluctance.  Already  the  same  year 
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Slater  '  made  a  careful  analysis  of  the  errors  that  could  be  expected  and 
showed  that  the  amount  of  agreement  obtained  with  experiment  was  by  no 
means  fortuitous.  Despite  that  the  Hartree  equation  was  not  generally 
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accepted  and  used  until  a  few  years  later  when  Slater  '  and  Fock  '  independent¬ 
ly  showed  that  the  Hartree  equation  augmented  with  an  exchange  term  could 
be  derived  from  the  variational  principle  using  a  Slater  determinant  as  trial 
function. 

The  HF  method  has  had  great  success  both  in  treating  atoms,  molecules 
and  the  band  structure  of  solids.  Concerning  the  band  structure  of  solids,  it 
is  really  not  at  all  obvious  that  the  HF  method,  considered  to  be  founded  on 
the  variational  principle,  should  give  sensible  results.  The  band  structure 
describes  the  excitation  energies  which  form  a  continuum.  An  excitation  energy 
is  obtained  by  subtracting  one  total  energy  from  another.  The  variational 
principle  tells  us  that  each  total  energy  is  the  best  possible,  which  means  it 
has  a  finite  and  quite  large  error.  Thus  the  variational  principle  furnishes  no 
direct  reason  why  the  HF  method  should  give  a  sensible  band  structure.  This 
kind  of  shortcoming  of  the  variational  principle  also  occurs,  to  a  lesser  extent, 
for  the  cohesive  energies  of  solids  and  for  the  excitation  energies  of  atoms 
and  molecules.  It  is  the  purpose  of  this  paper  to  try  to  explain,  from  a  many- 
-particle  standpoint,  why  and  when  the  HF  equations  give  a  reasonable  descrip¬ 
tion  of  excitation  spectra  and  also  to  describe  different  attempts  to  improve 
the  HF  equation  by  a  "correlation  potential"  or  '’polarization  potential". 

In  the  HF  method  the  total  wavefunction  is  described  by  one  Slater 
determinant.  A  Slater  determinant  is  really  a  very  poor  approximation  to  the 
true  wave  function  and  it  can  indeed  by  shown  that  the  overlap  integral  be¬ 
tween  any  Slater  determinant  and  the  true  wave  function  tends  to  zero 
exponentially  with  the  number  of  particles  in  any  normal  system.  In  a  Slater 
determinant  there  is  some  correlation  between  particles  of  the  same  spin 
brought  about  by  the  Pauli  principle  (the  Fermi  hole),  but  correlation  between 
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particles  of  opposite  spin  is  totally  missing* ' This  effect  can  be  studied  for 
e.g.  helium  where  very  accurate  wave  functions  have  been  calculated.  If 
one  keeps  the  position  of  one  electron  fixed  and  plots  the  wave  function 
against  the  distance  r  between  the  two  electrons  there  results  the  so 
called  correlation  hole  or  Coulomb  hole,  which  for  helium  is  a  curve  fair¬ 
ly  well  approximated  by  the  expression  C(l  +  r)  up  to  r  of  the  order  of 
one  atomic  unit. 

Though  this  correlation  feature  cannot  be  incorporated  in  any  Slater 
determinant  it  can  be  brought  into  the  HF  equation  itself  by  adding  a  correla¬ 
tion  term*  This  cannot  be  justified  if  one  looks  upon  the  HF  eq.  as  some¬ 
thing  providing  orbitals  for  a  complete  N  electron  wave  function  but  it  is 
very  reasonable  if  one  takes  the  original  intuitive  standpoint  Hartree  had 
when  first  deriving  the  equations.  Then  one  focuses  attention  on  one  specific 
electron  and  that  electron  should  not  see  the  average  potential  but  there 
should  be  a  correlation  hole  in  the  charge  distribution  of  the  other  electrons. 
There  are  several  treatments  of  such  polarization  effects.  The  simplest 
example  is  the  case  of  an  atom  where  one  electron  has  an  orbit  lying  outside 
the  remaining  ion.  The  polarization  of  the  ion  causes  an  extra  term  in  the 
potential  equal  to  -^al/r  ,  where  a  is  the  ion  polarizability.  Another 
famous  example  is  the  case  of  an  electron  moving  in  a  highly  excited  state 
round  a  charged  impurity  in  a  semiconductor.  The  energy  levels  are  then 
described  quite  accurately  by  a  hydrogen -like  formula  with  an  effective  mass 
and  a  dielectric  constant.  These  two  cases  have  been  justified  rigorously, 
the  first  by  van  Vleck  and  Whitelaw  and  independently  Mayer  and  Mayer 
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in  1933  and  the  second  by  Kohn  ’  in  1957.  The  single  particle  equations  in 
these  cases  are  really  equations  for  quasiparticles  or  dressed  particles, 
that  is  equations  for  a  bare  electron  and  a  polarization  cloud. 

It  was  more  difficult  however  to  justify  a  polarization  term  in  band 

8) 

calculations  of  crystals.  Attempts  were  made  by  Callaway  ^  in  1957  to 
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include  polarization  of  the  ion  cores  in  alkali  metals  and  by  Wohlfarth  '  in 
1950  to  use  a  screened  Coulomb  potential  in  the  exchange  term  of  an 
electron  gas.  It  is  only  quite  recently  that  one-electron  equations  with  a 
polarization  term  for  solids  have  obtained  a  proper  interpretation. 
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THE  BOHM-PINES  COLLECTIVE  APPROACH 

A  large  step  forward  towards  the  understanding  of  the  behavior  of 
electrons  in  a  metal  was  made  with  the  work  by  Bohm  and  Pines  in  1953, 
which  was  extended  and  refined  mainly  by  Pines.  The  basic  idea  was  to  treat 
first  the  collective  behavior,  i.e.  the  plasma  vibrations  of  the  electrons. 

The  plasma  vibrations  are  brought  about  by  the  long  range  of  the  coulomb 
force  and  are  responsible  for  the  long  range  correlation  between  the  elec¬ 
trons.  It  is  a  high  frequency  phenomenon  corresponding  to  energies  of 
typically  10  to  ZO  eV,  and  it  can  be  observed  by  measuring  energy  loss  of 
fast  electrons  passing  through  thin  metal  foils. 

As  coordinates  for  the  collective  plasma  vibrations  Bohm  and  Pines 

used  those  Fourier  components  of  the  electron  density  for  which  k  was 

smaller  than  some  cut  off  value  k^  .  They  also  kept  all  electron  coordinates 

and  thus  they  had  to  introduce  subsidiary  conditions  on  the  wave  function  to 

be  justified  in  treating  all  considered  coordinates  independently.  These 

subsidiary  conditions  have  caused  a  good  deal  of  trouble  and  form  a  weak 

point  in  the  Bohm-Pines  theory  as  also  does  the  uncertainty  in  the  choice 

of  the  cut  off  value  k  . 

c 

After  writing  the  Hamiltonian  in  a  form  exhibiting  the  collective 
coordinates  and  adding  terms  containing  collective  momenta,  the  effect  of 
which  was  nullified  by  the  subsidiary  conditions,  Bohm  and  Pines  performed 
a  series  of  canonical  transformations  which  brought  the  extended  Hamiltonian 
on  a  form  having  a  collective  part,  a  quasi-particle  part  and  supposedly  weak 
interaction  terms.  The  collective  part  had  essentially  an  harmonic  oscillator 
form’ and  the  quasi -particle  part  looked  like  the  original  Hamiltonicin  but 
with  the  coulomb  potential  replaced  by  a  screened  interaction  of  almost  the 
Yukawa  form,  e 

A  reasonable  first  approximation  for  this  Hamiltonian  is  a  product  of 
harmonic  oscillator  functions  of  the  collective  coordinates  and  a  Slater 
determinant  of  plane  waves.  The  collective  part  contains  the  long  range  cor¬ 
relation  and  the  Slater  determinant  can  be  improved  by  a  perturbation  ex¬ 
pansion  in  the  short  ranged  interaction  to  include  short  range  correlations. 
This  perturbation  expansion  is  conve  gent  contrary  to  an  expansion  in  the 
unscreened  coulomb  interaction,  1/  r  . 

In  calculating  the  total  energy  the  collective  part  is  essential  and 
brings  about  a  very  reasonable  result,  but  in  treating  the  single  particle 
behavior  the  plasma  part  can  be  left  out  in  the  first  approximation  since  the 
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large  excitation  energy  of  a  ‘^plasmon^*  makes  it  essentially  frozen  in, 
when  only  the  low  excitations  close  to  the  Fermi  surface  are  considered. 
The  quasi  particle  part  in  the  Bohm  Pines  Hamiltonian  has  been  success¬ 
fully  used  in  evaluating  different  single  particle  properties.  A  good 
review  of  the  earlier  calculations  is  given  by  Pines  in  Solid  State  Physics, 
vol.  1  (1955).  More  extended  calculations  have  been  made  by  Fletcher 
and  Larsson  (1958). 


THE  LINKED  CLUSTER  INFINITE  PERTURBATION  EXPANSION 

The  basic  physical  ideas  behind  the  single  particle  behavior  of  elec¬ 
trons  in  solids  is  well  described  by  the  Bohm-Pines  theory  but  more 
rigorous  developments  furnishing  both  exact  formal  results  and  quantitative 
estimates  have  utilized  field  theoretical  methods.  The  first  break  through 

for  infinite  order  perturbation  theory  treatment  of  many  body  systems  came 

12) 

with  the  work  by  Brueckner  '  in  1955.  Ordinary  Rayleigh -Schrbdinger 
perturbation  theory  gives  rise  to  terms  of  higher  and  higher  powers  in  the 
number  of  particles  N  ,  but  Brueckner  showed  explicitly  for  the  lowest 
orders  in  the  expansion  that  these  unphysical  terms  cancel  and  he  assumed 
that  this  mechanism  worked  for  all  orders.  That  this  really  is  the  case  was 

13) 

shown  by  Goldstone  '  who  proved  the  so  called  linked  cluster  theorem.  In 
the  Brueckner-Goldstone  expansion  every  term  is  represented  by  a  linked 
cluster  diagram  giving  a  contribution  at  most  proportional  to  N  .  Goldstone 
utilized  methods  and  results  from  quantum  field  theory  and  his  diagrams 
were  those  introduced  by  Feynman  in  1949- 

14) 

These  results  were  almost  immediately  followed  up  by  Hubbard  * 
who  published  two  papers  in  1957  on  collective  motions  in  terms  of  many- 
-body  perturbation  theory.  The  first  paper  was  purely  formal  and  has  become 
a  corner  stone  in  this  field  and  the  second  contained  a  numerical  calculation 
of  the  energy  of  an  electron  gas.  The  results  agreed  very  well  with  those  of 
Bohm  and  Pines. 
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THE  van  HOVE-HUGENHOLTZ  QUASI-PARTICLES 

Parallel  to  the  Brueckner-Goidstone-Hubbard  work,  van  Hove  and 
15) 

Hugenholtz  '  during  the  years  1955  to  1958  developed  their  own  technique 
to  handle  the  perturbation  theory  of  large  quantum  systems.  They  obtained 
results  which  in  many  respects  went  further  than  those  of  Goldstone  and 
Hubbard,  but  since  their  technique  was  much  more  complicated  it  has  not 
been  so  widely  used.  Their  work  centered  on  the  properties  of  the  low 
excitations  in  a  Fermi  gas,  and  the  most  important  results  are  perhaps  the 
following.  Consider  a  low  excited  eigenstate  of  an  Hamiltonian  with  no 
interparticle  interactions.  The  statevector  is  a  Slater  determinant  with 
N  particles  filling  a  Fermi  sphere  with  radius  kp  and  one  additional 
particle  close  to  the  Fermi  surface.  Switch  on  the  interactions  slowly 
(adiabatically).  The  resulting  state  is  then  an  asymptotically  stationary  state 
of  the  full  Hamiltonian  with  a  decay  time  proportional  to  l/(k-kp)  .  They 
also  showed  that  the  energy  of  the  perturbed  state  when  the  extra  particle 
is  on  the  Fermi  surface  is  equal  to  the  ground  state  energy  of  the  N  particle 
system  plus  the  average  energy  per  particle. 

The  van  Hove -Hugenholtz  results  apply  only  to  a  translationally 
invariant  system,  such  as  an  electron  gas,  and  the  assumption  that  the 
infinite  order  perturbation  expansion  really  converges  has  later  been 
criticized  by  van  Hove  himself.  Even  if  the  series  should  diverge  math¬ 
ematically  it  could  still  provide  a  good  asymptotic  representation.  The  re¬ 
sults  must  be  considered  as  a  large  step  towards  the  understanding  of  quasi- 
-particles.  Their  results  indicate  that  it  is  justified  to  use  a  model  for  the 
actual  system  with  a  single  particle  Hamiltonian  whose  eigenstates  are  in 
one  to  one  correspondence  with  slowly  decaying  states  of  the  true  Hamil¬ 
tonian.  Approximate  calculations  of  the  decay  time  give  a  mean  free  path 

of  the  order  of  100  000  A  at  room  temperature  for  an  energy  of  the  extra 
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particle  equal  to  kT  That  is  usually  much  longer  than  the  mean  free 

path  due  to  phonon  collisions  and  therefore  negligible.  The  results  are, 
however,  incomplete  in  the  respect  that  they  apply  only  to  the  case  with  one 
particle  outside  a  filled  Fermi  sphere  and  it  remains  to  show  that  a  model 
state  corresponding  to  many  excited  particles  has  a  true  energy  which  is 
the  sum  of  the  corresponding  quasi -particle  energies.  One  has  further  to 
show  that  the  adiabatic  transforms  of  the  low  excited  model  states  with  all 
possible  combinations  of  particles  and  holes  (unfilled  states  inside  the 
Fermi  surface)  really  furnish  all  low  excited  true  states. 


THE  LUTTINGER  TPIERMODYNAMIC  APPROACH 


A  more  general  approach  to  the  concept  of  a  Fermi  surface  and 

18) 

properties  connected  with  it  has  been  made  by  Luttinger  ^  in  a  series  of 
papers  published  in  I960  and  1961.  His  theory  is  applicable  to  non-zero 
temperatures.  The  partition  function  is  represented  by  an  infinite  aeries 

19) 

expansion  in  a  way  first  developed  by  Bloch  and  de  Dominicis  '  in  1958. 
The  expansion  is  formally  very  sin;ilar  to  the  B rueckner-Goldstone 
expansion,  the  Feynman  propagators  being  modified  to  include  temperature 
and  a  slightly  larger  class  of  diagrams  to  be  considered.  In  the  first  two 
papers  which  were  worke^d  out  together  with  Kohn  and  Ward  it  was  shown 
that  the  extra  diagrams  in  the  statistical  mechanical  expansion,  so  called 
anomalous  diagrams,  in  certain  cases  gave  rise  to  non-zero  contributions 
in  the  limit  when  T  goes  to  zero.  In  these  cases,  which  actually  include 
the  most  interesting  situations  such  as  a  periodic  lattice,  the  Brueckner- 
-Goldstone  series  has  to  be  modified.  During  the  detailed  investigation 
necessary  to  establish  this,  several  v^ery  interesting  results  appeared  as 
by-products.  Most  of  these  were  purely  mathematical  but  for  one  thing  they 
established  that  a  certain  class  of  quasi  -  stationary  states  of  the  interacting 
system  in  the  translationally  invariant  case  had  a  life -time  proportional  to 
l/(k-kp)  •  These  states  are  not  the  same  as  those  considered  by  Hugen- 
holtz  and  van  Hove  but  a-re  instead  related  to  the  Green  functions  or  single 
particle  propagators  used  by  Martin  and  Schwinger  and  others. 

In  three  papers  Duttinger  derived  some  very  interesting  results 
subject  only  to  the  restriction  that  the  infinite  order  perturbation  theory 
should  be  convergent  or  at  least  furnish  a  good  asymptotic  representation. 
He  found  that  n^^  ,  the  mean  occupation  number  for  the  plane  wave  state 
k  ,  had  a  discontinuity  on  one  or  several  surfaces  in  k- space  which  collec¬ 
tion  of  surfaces  he  defined  as  the  Fermi  surface  of  the  interacting  system 
In  the  translationally  invariant  case  the  Fermi  surface  is  determined  by  the 
relation 


(1) 


where 


is  the  chemical  potential,  is  the  enerry  in  the  non-interacting 

2  ^ 

k  /2  ,  and  Kj^(E)  is  the  real  part  of  the  so  called  self  energy  of  the 


case, 
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qua si -particle, 

Kk(E)  ±  i  Jk(E)  (2) 

where  we  have  plus  sign  when  E  <  p  and  minus  sign  when  E  >  p  , 

Jj^(E)  ^0  ,  Jj^(p)  =  0  .  The  energy  of  the  above-mentioned  quasi -stationary 
state  is 

Eh  £vc  (3) 

where  is  the  ground  state  energy  of  an  N -particle  system  and  Ej^  is 

determined  by  the  equation 

Ek  ^  Ku  (E|cj 

Jk(Ek)  is  the  inverse  of  the  lifetime  of  the  state. 

Luttinger  also  obtained  exact  results  for  the  low  temperature  heat 
capacity,  the  spin  paramagnetism,  the  compressibility  and  the  characteristics 
of  the  de  Haas -van  Alphen  oscillations.  Most  of  the  development  is  concerned 
with  the  translationally  invariant  situation  but  the  main  features  of  the  theory 
in  the  band  case  are  given  and  also  an  explicit  expression  for  the  heat 
capacity.  The  forms  of  these  expressions  are  very  similar  to  those  in  the 
non-interacting  case.  So  is  e.g.  the  heat  capacity  in  the  interactionless  case 
proportional  to  the  level  density  at  the  Fermi  surface  calculated  from  the 
orbital  energies  €  ^  whereas  in  Luttinger's  expression  the  quasi-particle 
energies  Ek  appear  and  the  integration  is  taken  over  the  true  Fermi 
surface.  Luttinger  also  showed  that  the  volume  in  k-space  inside  the  Fermi 
surface  is  the  same  in  the  non-interacting  and  interacting  cases.  This  means  , 
of  course,  that  the  Fermi  surface  is  unchanged  in  the  translationally  in¬ 
variant  case^  since  by  symmetry  it  has  to  be  spherical* 
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THE  GREEI'^  FUNCTION  FORMALISM 


Besides  the  pei’turbation-theoretic  treatment  of  many-body  systems 
there  is  a  large  branch  of  research  utilizing  Green  functions.  The 
Schrodinger  equation  is  then  converted  into  a  chain  of  coupled  differential 
equations  involving  higher  and  higher  orders  of  Green  functions.  The 
method  used  to  handle  this  system  is  usually  to  insert  a  specific  approx¬ 
imation  for  say  the  two-particle  Green  function  and  thus  closing  the  system 
of  equations.  Green  functions  or  particle  propagators  have  long  been  used 

in  field  theory.  Their  statistical -mechanical  versions  were  introduced  by 
2  i)  * 

Matsubara  '  in  1955  and  the  zero  temperature  forms  were  extensively 
treated  with  respect  to  their  spectral  properties  by  Galitskii  and  Migdal 
in  1958.  The  perhaps  most  well-known  paper  in  this  area  is  that  by  Martin 
and  Schwinger 


22) 


in  1959. 


We  will  not  review  any  specific  work  on  Green  functions  but  instead 

24) 

make  a  development  that  borrows  features  from  many  different  papers  ' 
and  leads  straight  to  the  derivation  of  a  single  particle  equation  where  the 
orbitals  and  orbital  energies  have  an  exact  interpretation  with  respect  to 
the  many  body  system.  The  orbital  energies  are  actually  identical  to  the 
quasi-particle  energies  used  by  Luttinger  and  the  orbitals  are  simply  related 
to  the  first  order  density  matrix  of  the  system.  Such  an  equation  was  first 
derived  by  Pratt  in  I960  starting  from  a  theorem  derived  by  Hugenholtz 
and  by  Bloch  and  Horowitz  Phillips  derived  a  similar  equation  in 

1961  from  Hubbards  theory  and  gave  a  discussion  of  the  order  of  magnitude 
of  the  different  terms. 


SPECTRAL  RESOLUTION  OF  THE  ONE-PARTICLE  GREEN  FUNCTION 
The  one -particle  Green  function  is  defined  as 

G  <N  1 T  ( V(»)  yY)  )  I  N  >  (5) 

where  x  stands  for  three  space  coordinates  and  the  spin  and  time  coordinates. 
i)j(x)  is  the  field  operator  in  the  Heisenberg  representation  and  |N>  is  the 
state  vector  for  the  ground  state  of  an  N -electron  system.  T  is  the  time- 
ordering  operator  placing  the  operator  with  the  latest  time  to  the  left  and  in- 
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serting  a  minus  sign  if  the  operators  and  change  places.  For  equal 
times  G  is  defined  as 

i  <Hl  *  i  W 

where  x  stands  for  space  and  spin  coordinates  and  p(x,  x’)  is  the  usual 
first  order  density  matrix.  G  has  a  very  direct  physical  interpretation 
which  becomes  clear  if  we  take  say  t  >  t'  and  rewrite  G  slightly  is 

the  ground  state  energy  of  the  N -particle  system)| 

5fe,»  is>')  -  -  kn  I  fy v>  I”  > '  -  i  N  ( rty  w  > 

H'.H-En  ,  T-  t-  .  (V) 

Thus  G  gives  the  probability  amplitude  that  a  particle  created  at  the  point 
X*  at  time  t^  reaches  the  point  x  at  time  t  .  That  statement  is  actually 
not  quite  true  since  the  states  involved  are  not  normalized  to  unity  but  the 
relation  anyhow  makes  it  reasonable  to  call  G  a  single  particle  propagator. 

We  will  now  consider  the  spectral  resolution  of  G(x,  x*;  t)  ,  where 
T  =  t  -  t*  , 


^  >o 


r  <o 


(8) 


where 


6ic  *  Eh+1,4  -Eh 


6^  «  Bh”  -C 


(9) 


In  a  HF  approximation  the  state  vectors  are  single  Slater  determinants  and 
the  g -functions  then  become  the  occupied  orbitals,  the  f -functions  the  virtual 
orbitals  and  the  €  -values  the  corresponding  orbital  energies.  In  general  there 
are  an  infinite  number  of  g -functions.  The  f -functions  and  the  g -functions 
are  further  not  in  general  normalized  to  unity  and  each  set  has  linear  depend¬ 
ences  between  its  functions.  We  now  make  a  Fourier  transform  with  respect 
to  time. 


-flo  k 


k 


(10) 
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To  have  a  unique  definition  for  all  c  an  infinitesimal  positive  convergence 
factor  A  has  been  introduced*  For  a  discrete  spectrum  G  has  simple 
poles  at  each  of  the  ’’quasi -particle"  energies  and  €  ^  .  For  a  con¬ 
tinuous  spectrum  the  sums  are  converted  into  integrals  and  the  poles  dis¬ 
appear. 


THE  SELF -ENERGY  OPERATOR  M  AND  THE 
QUASI -PARTICLES 

It  is  convenient  to  introduce  the  self-energy  operator  M(x,x’)  which 
can  be  considered  as  defined  by  the  relation  (the  Dyson  equation), 

G  =  +  G  MG  (11) 

o  o  '  ' 

where  G^  is  the  Green  function  when  there  are  no  interparticle  interactions 
and  G^MG  is  a  matrix  product,  x  and  x’  being  matrix  indices.  This  rela¬ 
tion  becomes  very  simple  when  there  is  translational  invariance  since  a 
Fourier  transform  with  respect  to  space  and  time  then  factorizes  the  matrix 
product  (G^  (^»  ®  )  =  *  "  all  matrices  are  diagonal  in  spin), 

Krj  Jj^»  ^  and  p  are  the  same  quantities  as  used  by  Luttinger-  M(k,  e  ) 
can  be  analytically  continued  to  complex  values  of  c  ,  the  continuations  being 
different,  however,  starting  from  the  part  of  the  real  axis  where  €  >  p  and 

the  part  where  €  <  p  .  We  now  consider  a  value  of  €  slightly  larger  than 
p  *  then  larger  than  zero  but  small.  We  may  then  assume  that 

G(k, €  )  has  a  pole  slightly  under  the  real  axis,  say  at  €  =  €  .  By  deform¬ 

ing  the  contour  of  integration  for  €  in  the  expression 

^  Me  ^ 

to  run  as 
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_ _ _ ^ 

_  .  _  v..  .  .  , 

\ _ 

■■  1 

V 

-  _  _  _ _ _ 

? 

we  obtain 

7  ^ - 7 - 

Fiji 


« *  {^6y,  *t  Kk  *" 

It  ^ 


(14) 


The  +  .  .  .  .  comes  from  the  contribution  for  e  =  [o,  +  i  e  ^  since  G(k,  e  ) 
is  different  in  the  upward  and  downward  paths,  these  paths  belonging  to  dif¬ 
ferent  analytical  continuations.  There  could  also  possibly  be  more  poles 

51? 

than  the  one  at  €  and  the  contour  has  then  to  be  deformed  correspondingly# 
Assuming  the  +  .  .  .  .  terms  to  be  small  '  we  find  that  G(k,  t)  describes 
a  quasi -particle  with  energy  ^  also  study  G  from 

its  spectral  resolution,  eq.  (8).  An  approximately  exponential  dependence 
of  T  will  result  if  the  k:th  Fourier  component  of  the  ff  -products  is  essen¬ 
tially  different  from  zero  only  within  a  narrow  energy  region  for  the  €  ^  . 

The  center  of  this  region  is  then  «  +  Kj^  and  its  spread  Jj^  « 


THE  EQUATION  FOR  G  AND  FOR  THE  AMPLITUDES 

V/e  now  return  to  the  case  with  no  translational  invariance.  After 
Fourier  transforming  the  Dyson  equation,  (11),  with  respect  to  time,  where 
we  still  have  translational  invariance,  and  utilizing  the  relation  G^ 

=  €  -  h  ,  where  h  is  the  relevant  single  particle  operator,  we  obtain 

(e  -  -fv  -  ^  1  (15) 

where  the  matrix  indices  x  and  x*  are  suppressed.  Inserting  the  spectral 
resolution  for  G(€  )  and  the  identity 

I  ^  ^ 

which  is  trivially  obtained  from  the  anticommutation  relation  of  i|j  and  ili^ 
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\ve  obtain 


-iLlLr\ 


tkC 


Now  let  €  approach  €  ^  and  assume  that  those  for  which  the  €  ^ 
are  equal,  are  linearly  independent.  We  then  have 


(e  ks.  (^H,) ^  f  vc  **  o  , 


This  is  the  single  particle  equation  we  set  out  to  derive*  The  assuinp.ivjn 
made  about  linear  independence  of  the  is  really  very  weak  as  can  be 
checked  by  considering  explicit  examples.  The  for  different  ^  ^  are 
on  the  other  hand  by  no  means  in  general  linearly  independent.  A  similar 
equation  is  of  course  also  obtained  for  the  .  The  equation  is  valid  when 
En+T  ^  ^  discrete  energy  level.  It  is  also  valid  when  ,  is  cont':  • 

uous  if  E  becomes  discrete  when  we  consider  only  a  specific  symmetry  of 
the  f  as  in  the  case  of  an  extra  particle  in  a  semi-conductor,  first  non 

trivial  approximation  to  this  equation  in  case  of  a  semiconductor  was  derivti: 
by  Pratt  '  (I960).  The  equation  in  its  general  form  was  suggested  by  Hedin 

and  Lundqvist  {i960)  but  the  amplitudes  g,  were  there  erronously  identi - 

29)  ^ 

fied  with  the  L6wdin  '  natural  spin -orbitals. 

For  a  discrete  spectrum  G(  e  )  and  G^(t  )  are  Hermitian  axid  thus 
M(e)  =  G^*^(c)  -  G  ^(c)  is  Hermitian,  but  for  a  continuous  spectrum  M(€) 
is  not  Hermitian  in  general.  We  can  then  consider  the  equation 


and  try  to  interpret  it.  When  the  M -operator  is  non -Hermitian,  the  eigv^n- 

values  in  general  become  complex.  In  analogy  with  the  analysis  made  for 

G(t)  in  the  translationally  invariant  case  one  might  expect  that  cp  correspond 

to  a  wave -packet  of  functions  and  the  real  part  of  t  to  their  mean  -^nargy 

In  case  of  one  extra  electron  in  an  insulator,  M  remains  Hermitian  for  a 

certain  range  of  €  -values.  This  case  has  been  treated  extensively  by  Kolin 
30) 

and  Ambegaokar  '  from  a  more  general  point  of  view  than  the  abo’-'e  equa¬ 
tion. 
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APPROXIMATIONS  TO  THE  M-OPERATOR 


To  make  the  equation  useful  we  now  have  to  obtain  approximations  to 

the  M-operator.  One  method  is  to  make  a  perturbation  expansion  for  M  , 

which  can  be  represented  by  a  well  known  set  of  diagrams#  We  can  also 

Z3) 

utilize  the  Martin  and  Schwinger  '  method  with  functional  derivatives,  and 
thus  circumvent  the  usual  convergence  difficulties*  It  is  then  possible  to 
obtain  an  expansion  in  higher  and  higher  orders  of  correlation  functions 


where  . 

^'(1)  =.  f[l)f(l)  -  <  ,  (2l) 

Each  order  is  obtained  by  an  iteration  procedure  and  the  expansion  could 
reasonably  be  expected  to  have  much  better  convergence  properties  than  the 
usual  one  in  •  The  details  of  the  derivation  are  evaluated  in  an  appendix 
since  the  possibility  of  using  an  iteration  scheme  is  not  pointed  out  earlier 
in  the  literature  to  the  knowledge  of  the  author# 

The  first  approximation  is  Hartree-like 


where  ^ 

V(s;  ~  <  fdi')  >  in'  ^  r{hi')  =  , 

and  <  P  is  the  true  density  of  the  N  electron  system.  The  next  itera 

tion  is, 

(file  -  -  V(ij)f<c(4)  '  /  6  k)  (24) 

where  od 

and  ^'l/(x,x*)  is  a  screened  potential  first  introduced  by  Hubbard, 

'\r  (jc,  %')  =  ^  x'j  h*)  (26) 

If  the  correlation  term  in  r  is  neglected  we  are  back  to  a  HF-like  approx¬ 
imation#  The  correlation  term  is,  however,  quite  tractable  since  it  is  re- 
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lated  in  a  simple  way  to  the  linear  response  function  of  the  system,  which 
in  turn  can  be  approximately  evaluated  from  a  time -dependent  Hartree  or 
HF  treatment.  G  could  be  approximated  with  the  solution  for  the  Hartree 
or  HF-like  equation.  The  next  term  is  much  more  complicated  and  hard  to 
deal  with,  though  Hubbard  has  tried  to  include  it  in  a  qualitative  way  in  his 
treatment  of  the  electron  gas-  Hubbard  treated  by  numerical 

methods  in  the  above  mentioned  approximation  (time  dependent  Hartree 
and  G  ==  G^)  whereas  Gell-Mann  and  Brueckner  (l957)  and  DuBois 
(1959)  made  an  expansion  in  r^  ,  the  interparticle  distance.  This  expan¬ 
sion  is  only  valid  at  most  up  to  r^  =  2  ,  but  one  could  hope  that  the  complete 
term  should  give  a  reasonable  approximation  also  in  the  metallic 
density  region  (r  =  2-5).  We  will  call  the  polarization  approxima¬ 

tion  for  M  ,  for  reasons  which  will  become  clear  in  the  next  section. 


DISCUSSION  OF  THE  POLARIZATION  APPROXIMATION  FOR  M 

The  formal  properties  of.  the  polarization  approximation  for  M  , 
eq.  (25),  such  as  its  spectral  resolution  and  its  relation  to  the  linear  response 
function  are  discussed  in  an  appendix.  In  this  section  we  will  utilize  these 
results  in  a  discussion  of  how  the  approximation  works  for  some  simple 
physical  systems  and  how  it  is  related  to  approximations  developed  by  other 
authors. 

Let  us  first  consider  the  simple  situation  of  an  electron  moving  well 
outside  a  closed  shell  ion.  The  energy  levels  are  then  to  a  first  approximation 
hydrogen -like.  Exchange  effects  have  little  influence  but  the  polarization  of 
the  ion  by  the  outer  electron  gives  a  small  shift  of  the  hydrogen -like  energy 
levels  which  has  been  observed  to  a  fair  accuracy  ^ 

A  simple  and  reasonably  accurate  way  to  calculate  this  shift  is 

34) 

furnished  by  an  "adiabatic"  approximation  i  The  outer  electron  is  then 
considered  as  a  fixed  point  charge  which  polarizes  the  ion.  The  induced 
charge  density  p(x)  gives  rise  to  an  extra  potential  <^(x) ,  which  acts  back 
on  the  outer  electron, 

jljL)  =  i:  ! ‘v'l&.’i'j  f(i‘)  (27) 

The  factor  ^  arises  since  -grad  <()(x)  has  to  give  the  force  on  the  outer  elec- 
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tron  and  the  gradient  acts  on  both  v  and  p(x')  ,  the  latter  containing  x 
implicitly.  Calculating  the  induced  charge  by  linear  response  we  have  from 
eq.  (Bl6) 


where 


r-'  fiOi)  Ps%') 


^  _  BN,0 

ps(^)  =  <N  , 


(28) 


(29) 


An  equivalent  result  can  be  obtained  by  a  quite  rigorous  treatment. 

Since  <(>(x)  is  a  small  perturbation  it  is  svifficient  to  consider  its 
expectation  value  with  respect  to  the  outer  electron  orbital.  Using  the  assump¬ 
tion  that  this  orbital  is  well  outside  the  ion  and  observing  that  the  expression 


z'- 


Bh,s  - 


(30) 


H,o 


is  essentially  different  from  zero  only  when  x  and  x'  are  more  or  less  in¬ 
side  the  ion  (the  expression  (30)  is  the  kernel  in  the  linear  response  function 
and  we  do  not  expect  much  induced  charge  outside  the  ion),  we  can  obtain  a 
good  approximation  of  <^(x)  by  making  a  multipole  expansion  of  the  coulomb 
potentials.  The  constant  terms  vanish  and  if  we  keep  only  the  dipole  terms 
we  have 


p-/  t  I  _  y* 

1 


JL  JL  1  Vi 

1^1'  Ix'i^  Z-t 


A: 


and  thus 


-  ± 

^  %  mi 


oC 


(31) 


(32) 


where  a  is  the  static  polarizability  of  the  ion.  This  result  can  also  be  ob 
tained  from  a  very  simple  argument  as  shown  by  Born  and  Heisenberg 
1924. 


in 


We  now  want  to  show  how  the  polarization  approximation  for  M  can 
reproduce  eq.  (28).  From  the  appendix,  eq.  (B4),  we  find  that  is  equal 

to  a  HF-like  exchange  term,  which  we  neglect  in  this  case,  plus  a  polariza¬ 
tion  potential  .  From  eq.  (BIO)  we  have, 

(f  k  ( ',^k.)  I fw)  -  (fk  I  T  f  k'  ft"  F(^k. f  k) 

It'  / 


(33) 
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where  is  the  outer  orbital  we  are  considering,  and  terms  with  \i 
have  been  neglected  since  they  contain  ion  orbitals  and  thus  are  small  where 
is  appreciable.  Since  the  ion  has  a  closed  shell  structure  we  can  assume 
that  the  ion  excitation  energies  are  much  larger  than  the  differences  €,  -  €,  * 
for  those  k'  where  fj^'  and  have  significant  overlap.  Thus  we  approxi¬ 
mate 


=  F(0)  IfO  = 


(34) 


and  have  by  that  obtained  the  same  approximation  as  the  "adiabatic”  one. 

Another  method  to  approximate  the  polarization  approximation  is  ob¬ 
tained  by  making  an  assumption  about  the  time  behavior  of  G  and  r  . 

y 

(t)  without  the  polarization  part  is  simply  v(t)  which  contains  a  6 -function 
in  T  .  We  now  assume  that  polarization  does  not  broaden  this  6 -function  peak 
too  much,  so  that  the  time  variation  of  G(t)  within  the  peak  still  can  be 
neglected  except  for  the  discontinuity  of  G(t)  at  t  =  0  .  Using  the  HF- 
-appi  oximation  for  G  we  have 


^cc  —1 

z)  =  (35) 

where  ^  TT  >  O 

^  I  o  T  <o 

Our  assumption  about  the  time  dependences  amounts  to  putting  e”^  ^  =  1  . 
Using  the  identity 

we  then  have 


From  (A25)  and  (37)  we  obtain 


(38) 
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Since  v(t^)  ©(t)  =  0  we  are  allowed  to  put  an  extra  v  into  the  first  term 
of  the  last  member  of  (38),  thus  emphasizing  that  this  term  only  contains 
the  polarization  part  of  y  .  Performing  a  Fourier  transform  with  respect 
to  time  we  have,  compare  with  (B6), 


0^ 


vr  -  L  (39) 


M(x*  x';  0)  thus  gives  the  "adiabatic*’  approximation  with  one  modification, 

namely  we  have  a  screened  potential  r  instead  of  v  in  the  exchange  term* 

This  modification  matkes  no  difference  in  the  ion  plus  outer  electron  case 

since  there  the  exchange  term  anyhow  is  small.  For  an  electron  gas  on  the 

other  hand  it  is  the  second  term  rather  than  the  first  which  is  important. 

The  first  term  is  local  in  space  and  thus  by  translational  invariance  it  is  a 

constant  and  gives  only  a  constant  shift  of  the  quasi -particle  energies.  The 

second  term  is  similar  to  the  one  used  by  Pines  to  evaluate  single -particle 

properties,  though  the  screening  lengths  are  somewhat  different.  Thus  e.g. 

the  anomaly  for  the  specific  heat  of  an  electron  gas  which  appears  in  a 

HF  approximation,  is  not  present  when  a  screened  potential  is  used,  as 

9) 

shown  already  by  Wohlfarth 


Although  the  first  term  is  trivial  for  an  electron  gas  it  may  be  quite 

important  for  a  real  solid.  This  term  has  been  considered  by  Callaway  for 

alkali  metals.  Unfortunately  for  this  discussion  he  only  estimates  its  effect 

on  the  cohesive  energy  rather  than  on  the  band  structure.  Phillips  '  argues 

from  a  different  type  of  approach  that  this  term  should  have  little  importance 

in  the  weakly  periodic  case.  It  should,  however,  be  realized  that  even  though 

this  term  might  be  essentially  constant  within  a  band  it  could  diffei- 

appreciably  between  bands.  Another  type  of  approach  using  a  "projected 

wave  field"  has  recently  been  made  by  Bassani,  Robinson,  Goodman  and 
37) 

Schrieffer  This  approach  is  adopted  to  utilize  the  rapid  convergence  of 

the  'Orthogonalized  Plane  Wave"  method.  The  basic  physical  idea  is  that 
the  core  electrons  can  be  treated  dynamically  independent  of  the  valence 
electrons. 

If  we  should  like  to  calculate  through  the  linear  response 

function  and  want  to  keep  the  approximation  (39),  we  have  to  take  it  in  the 
static  limit,  that  is  €  =  0  .  The  complete  expression  for  the  real  part  of 
Vp  is  given  by  (B24).  It  involves  the  real  part  of  the  response  function  as 
well  as  the  imaginary  part  which  is  a  kind  of  oscillator  strength, 
gives  a  measure  of  the  probability  for  a  transition  between  the  ground  state 
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-aX^uxs  ox  aouBADxaj  axx^TT  ^ABq  AxqBqojd  Aaqx  pun  tit  xJtBdjaxunoo  ou  aABq 

sx^XTqJo  xo^xa  asaqx  jo  auios  •  oo+  ox  oo-  uiojj  aSuBj  sanxBAuaSTa  ASjaua  aqx 
XBqx  puTj  aAv  lex's jado-x\  X^aptiadap- ABjiaua  aqx  qXT^^  uoTXBnba  aqx  jo  uinjcx^^ds 
aqx  japxsuoo  a/A  uaq^  •  joxHJado- jai  aqx  qXT^v  paxaauuoD  ‘poAxosun  jlvj  os 
•smaxqoad  x'BOxx'suiaqxBui  Aubui  ajB  ajaqx  X'sqx  paAjasqo  aq  pxnoqs  xi 


•ssajHojd  Jtapun  /aou  st  qjco/A  qang  ‘stnaxsAs  xaxdxxioo  ajoxu  ox  AxxxTR'^^TlddB 
aqx  ox  SB  aouBpinS  jaqxJnj  oatS  xTT^  xuaqx  uioaj  pauxBxqo  sxx^xsaj  aqx  ^sxjuox'b 
Xxsqx'E  aqx  puB  sbS  uoJxaax^  sb  sixiaxsAs  axduixs  uo  paxsax  aq  ox  uaqx  aABq 
SBapx  aqx  •suotx'Btxitxo JddB  x^s^xxoBJcd  x^o  qjoM  ox  Avoq  sb  apinS  Aub  UTBxqo  ox 
japjo  UT  ajnxaid  x'^aTsAqd  axduiTs  auios  qxT^^  suoTssajdxa  x'buijoj  aqx  xaauuoa 
ox  papaau  qDnui  AjcaA  st  xt  ‘uotxbtiijiojut  xtiBAaxaj  j:t  qanui  os  suTBXuoa  Jox'B;^ado 
-K  ^oj  uoTSsajdxa  xa^xa  aqx  aouTs  puB  aSjBX  os  st  luaxqojd  Apoq-AuBui 
aqx  JO  AxTxaxdTXioD  x'eoiX^u^aqxBUi  aqx  aauTS  x^qx  Xiaj  st  xi  •satxiaqas  uotx 
-Buitxo:rddB  poo8  paau  a/A  uotx'eituiijoj  '['SJauaS  aqx  puoAaq  paaaojd  ox 


•ASaaua  uoTXBXiaxa  ub  ppiA  ox  paxa'saxqirs  ajB  satSjaua 
TBXox  aSjB^  A'['[BaTdoDsoJOBTxi  o>ax  ajaqM  asBo  aqx  ut  paDUBTBqun  aaoui  st  uotx 
-BnxTs  aqx  SBajaq/A  saT^jaua  uoTXBXTDxa  aqx  ut  jLoaaa  tttsuis  AxSuTpuodsajaoa  b 
saATS  joxBjado  AS:taua  jjas  aqx  ui  jojcja  axBjapoiu  b  X'sqx  xaadxa  jcaqxJnj 
Astiz  a^  "UiaxsAs  Apoq-AuBixi  b  ut  os^b  SuiuBaui  x^Bxa  ub  sBq  qDiqM  uoTxdxjasap 
B  UT  ‘auiTX  ajT]  B  puB  ASaaua  ub  SuiABq  ‘apxxJBd-TSBnb  jo  ‘a^oix^^d  b  jo 
xdaouoD  -[BDTSSBXO  aqx  asn  ox  a^qB  ajB  a/A  snqx  'SBapx  xBOisAqd  a^duiis  qXT^^ 
uoixoauuoD  asojo  sxt  si  uoixounj  aAB/A  uojxaa^a- Aubui  b  Aq  uoixdTjosap  b  qxjAV 
pajBduzoa  sb  uisijbuzjoj  uoixounj  uaajr)  aqx  jo  aSBxuBApB  utbxu  aqx 


SMHVT^3H  ONiamONOO 


•(uoTXBOTuruxiuioD  axBATJd  'xX^^d  'M'O) 
uoTXBZTJBTod  auios  UTBXUOD  AxTBaj  pxnoqs  uijax  XT^^o-uxds  aqx  osxv  •JoxBjado 
Yi  aqx  ox  pappB  aq  A^duiTs  ubd  ajnxanjxs  jBnsn  aqx  jo  uijax  xiq^o-uids  b 
AjBssaaau  jj  ‘axBxs  punoj8  aqx  aAoqB  co  ASjaua  ub  SuiABq  ax^xs  paxtaxa  ub  pUB 
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h’.iF.  been  carried  out. 

Also  in  the  present  context  we  are  faced  with  a  self-consistency 
problem.  When  a  specific  type  of  approximation  for  the  M -operator  is 
chosen  we  have  to  iterate  to  obtain  self-consistently  the  orbitals  which 
form  the  propagator  G  .  In  a  specific  step  of  this  procedure,  solving  for 
the  orbitals  using  a  rough  potential,  we  may  of  course  utilize  the  varia¬ 
tional  principle.  We  may  also  convert  this  intermediate  problem  to  ma¬ 
trix  form,  expanding  the  orbitals  in  a  fixed  truncated  basis  set  of  say 
r  functions.  The  r  orbital  solutions  we  then  obtain  provide  the  corre¬ 
sponding  truncation  of  the  infinite  summation  appearing  in  ,  eq.  (B24). 

The  present  formalism  is  so  far  inadequate  in  treating  multiplicity 
sTjlittings.  The  main  applicability  is  to  atoms,  molecules  and  solids  with 
one  electron  outside  a  closed  configuration  and  to  metals  where  there  are 
JO  many  electrons  outside  a  closed  configuration  (the  ions)  that  coupling 
between  spins  of  different  conduction  electrons  become  unimportant- 

The  instability  sometimes  observed  in  calculating  orbital  energies 
in  the  HF  method  could  very  well  be  related  to  this  failure  of  taking 
multiplicity  splittings  into  account.  The  density  matrix  is,  however,  quite 
stable  and  that  is  all  that  is  needed  as  regards  the  excited  orbitals. 
Specifically  as  regards  band  structure  the  exact  theory  makes  predictions 
only  for  excitations  close  to  the  Fermi  surface  and  there  is  hence  an  addi¬ 
tional  reason  to  regard  the  lower  bands  as  less  significant. 

It  should  finally  be  recognized  that  the  development  given  in  the  last 
few  sections  of  one -electron  behaviour  is  by  no  means  complete.  Thus  e.g. 
the  vibrations  of  the  atoms  can  very  well  give  an  influence  on  the  quasi - 
-particle  spectrum  comparable  to  the  polarization  effects  Further  in 

considering  transport  processes  not  only  the  quasi-particle  spectrum  is 
important  but  also  the  effective  field  inside  the  solid,  the  calculation  of  which 
is  quite  a  problem  in  itself.  Transport  processes  are  perhaps  most  effective¬ 
ly  treated  with  the  temperature  dependent  Green  function  formalism. 
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APPENDIX 

A.  APPROXIMATION  OF  THE  M  OPERATOR 
THROUGH  AN  ITERATION  PROCEDURE 

The  equation  of  motion  for  the  field  operator  in  the  Heisenberg  repre¬ 
sentation  is 

■i, »  VW' [ yWjHl .  (Ai) 

After  evaluating  the  commutator  we  have 

(ih-i fe)  -  i h'J  ==.  o ^  (A2) 

where 

Ip) ^i^)j  we'" ^  fl-fc 02 

The  functions  \x^  form  a  complete  orthonormal  set  and  a^^  is  an  annihila¬ 
tion  operator.  An  equation  for  the  Green  function  is  now  easily  obtained  from 
(A2). 

where 

Gii^)  =  (-0 

^  <T  fy  (1)  H'(^)  ^Hl'jy 

(i"^)  ^  O  dnJ  infiniffSirna./. 

In  eq-  (A4)  the  first  and  second  order  Green  functions  are  coupled.  Martin 
23) 

and  Schwinger  '  have  shown  that  the  term  with  G(l,Z;  1\2')  can  be  ex¬ 
pressed  as  a  functional  derivative  of  G(l,2)  .  Since  their  proof  is  rather 
complicated  and  not  very  explicit  we  will  rederive  the  functional  derivative 
expression  by  elementary  methods. 


(A4) 


(A  5) 


Consider  a  Hamiltonian  H  in  the  Schrbdinger  picture 
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Hi  =  Hs 


(A6) 


where  H  is  identical  with  the  Hamiltonian  H  defined  in  (A3)  and  W  (t) 
s  s 

is  a  small  perturbing  time  dependent  potential 


W  =  J  f(x)  h  ^ 


(A7) 


where  w(x,  t)  is  taken  as  spin  independent  (x  involves  in  general  both  space 
and  spin).  The  time  evolution  of  the  state  vector  in  the  Schrtidinger  picture 
is  described  by  an  operator  V(t,t'), 

[t>j  =  |t'>s  ,  (AS) 

Since  satisfies  the  Schrbdinger  equation,  V(t,  t')  has  to  satisfy 

i  h  v(M'J  =  (nf  -  (h)  V(t,  i') )  v(i,  f j  =  1 . 


We  also  introduce  the  unperturbed  time  evolution  operator  U(t,t')  =  e 
which  satisfies 

The  differential  eq.  (A9)  is  then  converted  into  an  integral  eq. 


-iH(t-t') 


(A  10) 


(All) 


The  crucial  quantity  we  need  is  the  functional  derivative  of  V(t,t')  with 
respect  to  the  perturbing  potential  w(x,  t)  ,  which  from  (All)  becomes, 

SVM  _  SVd.i'}  ^ 

S-^D.)  -  f'w'ffc/,.)  ,  ,  . 

V  — t  (An) 

f*1  i>o 
t<c 

^  .  6V  . 

when  is  in  the  time  interval  determined  by  I  and  t*  ,  otherwise  is 

zero.  Putting  w  =  0  in  (A12)  we  have 

The  general  expression  for  the  functional  derivative  in  question  is 

^  a- 1')  y(t, J  y{ix,  tj  , 


(A  13) 


(A  14) 
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as  can  be  verified  by  substituting  (A14)  into  (A12). 

It  is  convenient  to  introduce  the  following  type  of  expectation  values 
to  use  in  the  Green  functions 

where  il/pj(x)  =  V(t^,  t)  i|/(x)  V(t,  t^)  •  We  consider  t^  to  be  some  time  in  the 
distant  past  and  to  be  the  ground  state  of  H  .  T  is  a  time  far  in  the 

future  and  we  do  not  allow  w(t)  to  be  different  from  zero  for  other  t  than 
those  in  the  interval  between  and  T  .  When  w(t)  is  zero  for  all  t  the 
definition  (A  15)  gives  the  expectation  value  of  Heisenberg  operators  with 
respect  to  the  ground  state,  cf  eq.  (5).  The  matrix  element  in  (A  15)  has  the 
same  structure  as  those  used  in  scattering  theory,  that  is  a  matrix  element 
between  unperturbed  states  at  times  far  in  the  future,  U(T,  ^ 

far  in  the  past,  1  '  with  a  scattering  matrix  V(T,  inserted. 


With  the  present  definition  of  matrix  elements  the  equation  for  G(l,  1’) 


becomes 


The  term  containing  the  second  order  Green  function  can  be  rewritten  by  aid 
of  the  equation 

-=  -  Cr(f(2.)  v>(i) 

,  (A17) 

=  -  +  i' f(%)y . 

Thus  equation  (A  16)  becomes 

(if*  -  lO)  -  V(i))<S.(V'j  - 1 


where 


V(ij  =  ^  W  +  J‘'r(i,?L)<f(2,)>  i(l),  (A19) 

To  obtain  an  iteration  procedure  for  solving  (A  18),  it  is  essential  to  introduce 
the  inverse  of  the  Green  function, 


(A20) 


J  (S.''’(1,%)S(5.,1')  Hi.)  =  ^(v‘)  ■ 

We  further  note  that  ^ 


i  'W'i'b} 


(A21) 
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We  choose  to  define  the  zero  order  Green  function  by  the  Hartree-like  equa¬ 


tion, 


(A22) 


(i  ft  -  W  -  V  (1) j  <? (1,1' )  =  f  (I, !’)  ^ 

and  hence  the  self  energy  operator  M  becomes 


(A25) 


since  (A23)  inserted  in  (A  18)  gives 

With  a  specific  approximation  for  M  ,  the  corresponding  approximation  for 
G  can  be  obtained  from  eq.  (A24)  and  put  back  into  eq.  (A23),  thus  giving 
the  next  iteration  step  in  the  approximation  sequence  for  M  .  The  first 
approximation  for  M  is  obtained  by  using  G^^^  in  eq.  (A23), 

=  -  i  /  7(1, H)  im 

=  ^  7(1.%} 

Here  the  Hubbard  effective  potential  Y  has  been  introduced, 

7(1, X)-  V(a)  m  = 

=  Tfi^)  -  i  !ir(o)  O' f'(i(j)y  Y(‘iix) 

The  next  approximation  for  M  is  obtained  by  calculating  from  the 

equation 

(if -  J fi^'\i.x}&%,i')  1(1)=^ SM,  (A«) 


(A26) 


and  then  use  the  solution  in  (A23), 


J^)  Gi'  ■  (i'2j  &') i&)  = 

=  i  &\i.}7(i,i}  +  i  ^  kk 

To  calculate  the  functional  derivative  of  we  need  those  of  G^^^  and  y. 


i 
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= -/v(4;v;<j,s;<T('/’i'};ff^;^f^yj>  . 

(2) 

The  expression  for  M'  '  becomes, 


(A29) 


If  we  use  the  symbols 


1  2*  -f  JL 

Yftij  -v-fti^ 


(A31) 


the  approximation  we  have  now  obtained  for  V+M  can  be  written  (we  neglect 

and 


the  difference  between  G 
for  G)» 


1.1 


1. 


and  represent  them  with  the  symbol 


(A  32) 


(3) 

where  p'  '  stands  for  the  three  particle  correlation  function  in  eq.  (A 30). 
The  first  term  in  and  obviously  converges  to  iG(l,  2)^V^(l,  2), 

and  we  have  not  always  been  careful  to  state  what  kind  of  approximation  for 
G  is  considered*  The  perturbing  potential  w  is  of  course  to  be  put  equal 
to  zero  once  the  expression  for  the  required  approximation  has  been  obtained* 
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B.  THE  SPECTRAL  RESOLUTION  OF  AND  ITS 

RELATION  TO  THE  LINEAR  RESPONSE  FUNCTION 


We  Fourier  transform  eq.  (A24)  with  respect  to  the  time  difference, 

t-t'  =  T  , 

[c-  (bi) 

For  the  M  operator  we  take  the  approximation  (A25) 

0D  ^ 

The  v(l^,2)  term  in  Y  (1,2)  ,  eq.  (A26)  gives  a  HF  type  exchange  term. 


For  the  remaining  part  of  the  operator  we  introduce  the  name  polariza¬ 

tion  potential,  V  (1,2;€)  , 

The  spectral  resolution  of  G  has  already  been  given,  eq.  (10).  We  rewrite 
it  in  a  slightly  more  compact  form 


G  6^  —  X. 


^  —  ^/T7 


where  cp  stands  for  both  the  f,  and  the  g-  and  €  is  either  €,  -  iA 
^  m  k  m  k 

or  +  iA  .  Since  €j^  >  p  >  the  poles  of  G(l,  2;  c)  are  distributed 
according  to  figure  2  . 

Figure  2. 


Note  that  ^  is  somewhat  arbitrary  for  a  discrete  spectrum  and  could  there 
be  taken  as  e.g.  the  mean  value  of  the  smallest  and  the  largest  . 
We  also  make  a  spectral  resolution  of  the  part  of  Y  that  remains  after  v 
is  taken  awayi 
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•• 

•• 

*_J  e  ^  T  {'i)  S  '^(M)  L  f(^f)y  ^(^j  ^ 

+  <  "" f>(&-»)y  0(-r:)  -  <  f(^-^  >  <  f^i)y']  == 

_  r'  f’i%,)fikz}  r^, 

L  -T777ZT2  ^  nts-}^)  niZl'^J, 

_- ,  r  <^  I /  /V.  5 > 


^6  -  ^N,S  -  JEn 


The  prime  on  the  sums  stands  for  deleting  the  term  with  s  =  0  .  In  writing 
the  last  step  of  (B6)  we  have  used  the  fact  that  ^  ^s(^J  is  real* 

which  is  true  also  when  we  use  complex  wave  functions.  Thus  we  have  for 
the  polarization  potential  V  f 


We  first  consider  the  evaluation  of  a  typical  term  of  the  contribution  from 
F(1,2;€')  , 


m 

I 


e-6' 


m 


u 


(B8) 


We  can  close  the  contour  by  a  Large  semi  circle  in  either  the  upper  or  the 


m-2 


We 


lower  half  plane  since  the  integral  tends  to  zero  as  fast  as  (c*) 

choose  the  upper  half  plane  in  this  case  and  thus  we  have  a  contribution  only 

when  €  corresponds  to  a  particle,  €  =  €,  • 

m  ^  ^  ’  m  k  ' 


«£ -<£m  -  <£i  i;4  0 


(B9) 


The  imaginary  part  of  can  now  be  absorbed  in  iA  and  we  consider 

^m  henceforth.  The  contribution  from  F(l,2;  -c')  is  obtained 


similarly  by  closing  the  contour  in  the  lower  half  plane.  The  complete  ex¬ 
pression  for  Vp  becomes, 

=  1 w A-'W  w  p/w 


~  ^  ‘fm  W  C  ^(lX)£- -£m)  Ofdm-/*)  -  ^(lX)  jJ  , 


(BIO) 
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We  now  want  to  relate  the  polarization  potential  V  to  the  linear  response 
function  R  .  To  obtain  R  we  consider  a  small  external  charge  density, 
p^^(x,  t)  ,  and  calculate  in  a  linear  approximation  the  induced  charge 
density  This  induced  charge  density  gives  rise  to  an  induced 

potential  <^(x>  t)  •  We  define  R  by  the  relation 


-  i;  I ijc',  t) U' 


The  induced  charge  density,  ,  from  a  perturbing  potential  w  ,  is 

according  to  (A8)  and  (All),  in  the  linear  approximation. 


f  I  f(^)  v(t,u)luy  -  <UluMf>(ij  - 

t 

Ool  d*' .  (B»2) 

We  take,  cf  eq.  (A7), 

f fit')  ftx f  (f,  0  k' k* .  ^ 

From  (B12)  we  then  have,  putting  t^  =  -oo  , 

f  f 

and  comparing  with  (B  1 1)  we  obtain 

R  (i.x)i,-h)^  -  i  J  <<,|  C  f(i'),  f(t:)]luy  viijjvltx) 

When  we  compare  R  with  the  polarization  part  of  y  ^  eq*  (A26)»  we  find 
that  the  main  difference  is  that  R  has  a  commutator  instead  of  a  time- 
-ordering  operator-  The  Fourier  transform  of  R  with  respect  to  time  is, 

R ^,1) -  I  ^ =  if -f- 

The  reason  that  we  have  taken  0  to  oo  as  integration  limits  rather  than  -oo 
to  00,  is  that  R(l,2;  t)  does  not  appear  for  t  <  0  in  eq.  (Bll).  We  now 
relate  the  imaginary  part  of  Vp  to  R  •  From  the  definition  of  F  we 
have  that, 

1.)*;  -  o 


for  €  <  0  , 


(B17) 
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and  thus  from  (BIO)  and  (Bl6), 


(B18) 


»pVi‘-;‘5/=‘  'p 
where  we  have  used  the  relations, 


^  l  for  e  >  0 


(B19) 


From  (B18)  we  find  immediately, 

^  O  ^  (B20) 

which  relation  has  to  be  satisfied  from  general  considerations,  as  shown  by 
Luttinger.  We  next  want  to  relate  the  real  part  of  to  R  ,  and  to  that 
end  we  write, 

F%V,i)  =  .  (B21) 


The  second  term  of  (B21)  can  be  expressed  in  R 


(2) 


Since 


we  have 


/  ^ 

s  ^  ^ 

From  (BIO),  (BZl),  and  (B23)  we  finally  have. 
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